Quantum Phase Transitions in Josephson-Junction Chains 
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We investigate the quantum phase transition in a one-dimensional chain of ultra-small supercon- 
ducting grains, considering both the self- and junction capacitances. At zero temperature, the system 
is transformed into a two-dimensional system of classical vortices, where the junction capacitance 
introduces anisotropy in the interaction between vortices. This leads to the superconductor-insulator 
transition of the Berezinskii-Kosterlitz-Thouless type, as the ratios of the Josephson coupling energy 
to the charging energies are varied. It is found that the junction capacitance plays a role similar 
to that of dissipation and tends to suppress quantum fluctuations; nevertheless the insulator region 
survives even for arbitrarily large values of the junction capacitance. 

PACS Numbers: 74.50.+r, 67.40.Db 



Quantum phase transitions, which are induced by 
quantum fluctuations at zero temperature, are distin- 
guished from classical phase transitions in several impor- 
tant respects; this has attracted much attention in recent 
years jjj. In particular, advances in fabrication tech- 
niques have made available arrays of ultra-small super- 
conducting grains, where the charging energy dominates 
the Josephson coupling energy and accordingly, quantum 
fluctuation effects are of paramount importance. Such 
Josephson-junction arrays have become a prototype sys- 
tem displaying quantum phase transitions between the 
superconducting and the insulating phases. In the vicin- 
ity of the superconductor-insulator transition, the fluc- 
tuation effects depend crucially on the dimensionality 
of the system. In the case of two-dimensional (2D) ar- 
rays, rich effects of quantum fluctuations and resulting 
phase transitions have been examined for a rather gen- 
eral form of the capacitance matrix, although there still 
exist unsettled issues in the quantum regime, such as low- 
temperature re-entrance On the other hand, one- 
dimensional (ID) chains of Josephson junctions, where 
quantum fluctuations should be more important, have 
been studied mainly in the two limiting cases: the self- 
charging model and the nearest-neighbor model where 
only nearest neighboring charges interact In the 
ID system with only self-capacitance, the role of dissipa- 
tion on the quantum phase transition as well as the 
persistence current and voltage || has also been consid- 
ered. 

This paper investigates the quantum phase transitions 
in general Josephson-junction chains with both the self- 
and junction capacitances. At zero temperature, the sys- 
tem is transformed into a two-dimensional system of clas- 
sical vortices, where the junction capacitance introduces 
anisotropy in the interaction between vortices. This 
leads to the superconductor-insulator transition of the 
Berezinskii-Kosterlitz-Thouless (BKT) type §,@, as the 
ratios of the Josephson coupling energy to the charging 
energies are varied. Interestingly, the junction capaci- 



tance here plays a role similar to that of dissipation and 
tends to suppress quantum fluctuations, enhancing su- 
perconductivity. However, the suppression is not strong 
enough, and the insulator region still remains even for 
arbitrarily large values of the junction capacitance. 

We consider a one-dimensional array of N ultra-small 
superconducting grains, each of which is coupled to 
nearest-neighboring grains via Josephson junctions. The 
system is characterized by three energy scales, Josephson 



coupling energy Ej and charging energies Eq 



72Co 



and Ei = e /2Ci, where Co is the self-capacitance of 
each grain and C± is the junction capacitance between 
nearest-neighboring grains. It is described by the Hamil- 
tonian 



H = 
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N 
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where the number n, of Cooper pairs and the phase 
(j>i of the superconducting order parameter at site i are 
quantum-mechanically conjugate variables: [ni,4>j] = 
iSij, and the energy has been rescaled in units of the 



Josephson plasma frequency Uui p 
have introduced Kq 



V8E Ej. Here we 



yj Ej 1 8Eo to describe the case 
Co =/= 0, and written the capacitance matrix in the form 



C;. 



A 2 (2Si 



<5i, 7 — 1 



where A 2 = Ci/Cq. On the other hand, when Co = 0, it 
is convenient to introduce Ki = d Ej / 8Ei . Note that A 
can also be written in the form: A = Ki/ Kq. 

In the imaginary-time path-integral representation, the 
partition function of the system reads 
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with the action 
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where the temperature has been rescaled according to 
hujpP — > /3. In Eq. (||), 9^ and d T denote the difference 
operator with respect to the position x and to the imag- 
inary time r, respectively: d x 4>j. T = 0j+i, r — 4>j,T arid 
d T (f>j yT = 4>j, T +\ — 4>j,T, arid the (imaginary) time slice St 
in the interval [0, (3] has been chosen to be unity (in units 
of huip) [|IlJ . Our main concern here is the quantum phase 
transition at zero temperature, and accordingly, the limit 
f3 — > oo as well as the thermodynamic limit N — > oo is to 
be taken. We then apply the Villain approximation fl2] | 
and integrate out {4>j, T }, to obtain the partition function 
in terms of an additional set of integer variables {wj jT } 
(as well as {nj, T })- The corresponding action is given by 



L[n, m] 
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where the two sets of integer variables satisfy the con- 



straint d r m 



3,T 



d T rij T 



0. This constraint is con- 



veniently taken into account by introducing an integer 
field A j - defined on the space-time dual lattice (j,f) = 
(j + 1/2, t + 1/2) in such a way that m J>x = —d T A^ - 

and ua T = d x A~, ~. Henceforth we will work on the dual 

j ' j 

lattice, and drop for simplicity the tilde sign over site 
indices. The partition function is thus written in terms 
of the unconstrained summation of exp{— over 
Aj tT 's, with the action 



L[A] = ^-J2( d * A ^)C- J 1 (d x A J , T ) + -±-^( dr A j>T ) 



2K 
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Finally, the Poisson summation formula, which decom- 
poses Aj jT into a real-valued field and a new integer field 
Vj. T , allows us to integrate out the real- valued field A^ T . 
This leads, apart from the spin-wave part, to the 2D sys- 
tem of classical vortices described by the Hamiltonian 



H v = 2ir 2 K Q ^ v i,t U (i-j, t~<j) v jiCn 
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where the vortex interaction is given by 



U(x,t) 
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with the Fourier transforms of the lattice Laplacian 
A(z) = 2(1 — cosz) and of the capacitance C(q) 
I + X 2 A(q). Note that unless Co = 0, the diagonal piece 
Uq = U (x=0, r=0) becomes arbitrarily large as N and 
/3 are increased to infinity. Therefore the diagonal term 
2ir 2 K Uo(J2j t v 3, t ) 2 snou ld vanish in Eq. (||), leaving 

j T = 0, and the 



the vorticity neutrality condition: y\ r - 
Hamiltonian (JfJ) may be written in the form 



IT ,J(T 



with U(x,t) = 2tt[U - U(x,r)]. The behavior of the 
(reduced) vortex interaction U(x,t) for several values of 



A is displayed in Figs. [I] and |2| which manifests the log- 
arithmic behavior at large length scales. In particular, 
the short-range anisotropy prominent for large values of 
A decreases rapidly with the distance. According to the 
renormalization group (RG) theory of the 2D Coulomb 
gas, such short-range anisotropy should not affect the 
universality class. Thus in the spirit of the RG theory, 
the system is expected to exhibit qualitatively the same 
critical behavior as the 2D Coulomb gas jl3) . 

We now investigate the quantum phase transitions dis- 
played by the Hamiltonian in Eq. (||). First we consider 
the simplest case of C\ = 0, which has been studied in 
Rcf. B. In this self-charging limit, the interaction in 
Eq. (0) takes the simple form U(q, us) = [A(g) + A(w)]~ 1 
or U(x, t) w In y 'x 2 + r 2 + | In 2+^e, where is the Eu- 
ler number. Thus the Hamiltonian (||) describes precisely 
the isotropic 2D Coulomb gas, and the system undergoes 
a BKT transition |9|, [l0| , |l2| from the insulating phase to 
the superconducting one as Kq is increased 0] . 

We next consider the opposite limit Co = 0, where the 
relevant (dimensionless) coupling constant is K\ instead 
of K . Accordingly, it is proper to use hoj p = ^8EiEj 
in rescaling the energy and the imaginary time, which in 
turn gives C(q) = \~ 2 +A(q) and [/(<?, w) = [1 + A(w)] -1 . 
As a result, we obtain the Hamiltonian in the form 



H v = 2ir 2 K 1 



lT,j<T 



(9) 
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where, in sharp contrast to the previous self-charging 
case, the vortex interaction U(x, r) = 2tt[Uq — U(x, r)] 
is short-ranged: 



U(x,t) 



2ir 
7! 



1 - S x ne' 



-\r\ 
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for |t| 1. It is further of particular importance that Uq 
in general does not diverge, which implies that the diag- 
onal term in the Hamiltonian (^|) does not give the vor- 
ticity neutrality condition. Thus (unbound) free vortices 
become pervasive, and the system remains insulating for 
any nonzero value of K\ . This can also be understood in 
the charge representation, where the partition function 
reads 



( 9> w)| a G- 1 (?,a;)L (11) 



with G(q, u>) = A(q)[A(w) + It follows from the an- 
alytic continuation of G{q,u>) that the charge excitation 
has a gap E g ~ huj p . Thus the long-range interaction 
between charges gives rise to the gap, resulting in 
an insulator. It is also of interest to compare this case 
with the "nearest-neighbor model" in Ref . || : In the lat- 
ter, only the nearest-neighboring charges interact, which 
allows only bound vortices. Consequently, the system 
is superconducting for all finite values of the interaction 
strength. 

We now turn to the more realistic case that A is small 
but nonzero, i.e., Cq 3> C\. In this case, the vortex 



interaction U is, to the order of A 2 , isotropic and has 
the asymptotic behavior: U(x,t) ~ In \J x 2 + t 2 + e(A), 
where e(A) is the vortex pair creation energy (per vortex) 
given by e(A) « (3/2) ln2+ 7i5 + (7r-l)A 2 + 0(A 4 ). Hence 
the system again reduces to the 2D Coulomb gas, with 
the vortex fugacity y(X) = exp[— 2irKoe(\)] diminished 
by the factor exp[— 2ttK (tt — 1)A 2 ]. The standard RG 
theory of the BKT transition jl0],[l2| then predicts the 
transition point K§(\) slightly decreased from K§(0) by 
the amount 



6K$_ 2^(0)-4 A2 
K§(0) y > 2ttKS(0) - 3 



0(A 4 ) 



T.6A 2 . 
(12) 



In this limit (A <C 1), the Hamiltonian in Eq. (0) can 
also be represented in terms of phase variables, which 
yields a 2D XY model with an additional interaction for 
nonzero A. To see this, we neglect C(A 4 ), and write the 
inverse capacitance matrix in the form 



(1-2A 2 )^+A 2 (<5, J+1 +^_ 1 ) 



(13) 



Note that without the off-diagonal term in , the first 
and the second terms in the action (|^) would just be 
the Villain form of the cosine action along the r direc- 
tion HH . To examine the effects of the off-diagonal term, 
we use the identity 



exp 



4K, 



= — = / dz 



exp 



which allows to separate the charge variables at dif- 
ferent sites. The resulting action then takes the Vil- 
lain form of the cosine action with the argument 



d T <t>j,T+i\J 'A 2 / Kq d x Zj,r along the r direction. To the 
order of A 2 , the Gaussian integration over Zj jT leads to 
the effective phase Hamiltonian 



H cS [<f\ = -Kq 



1 



A 



y^[cos(0a!0j,r) + COs(d T ^j, T )] 



X 2 K ^2[d T sind x 4> jtT ] 2 , 
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where the anisotropy has been removed by rescaling again 
the r axis by the factor (l-A 2 )- 1 /^ 2 /^. Note that 
Eq. (|H) reduces, for A = 0, to the standard 2D XY 
Hamiltonian. Here the junction capacitance not only en- 
hances the effective coupling of the XY model but also 



introduces an additional interaction given by the second 
term in Eq. (Q). Interestingly, the latter is very similar 
in form to the dissipation term in the effective action, 
which is known to suppress quantum fluctuations HQ. 
Therefore both effects contribute to the enhancement of 



3 



phase coherence, and it is concluded that the junction 
capacitance in general tends to suppress quantum fluctu- 
ations induced by the self-capacitance, thus reducing the 
insulator region. 

For larger values of A, Figs. |] and ||(b) show that the 



vortex creation energy e(X) and the short-range behavior 
of the interaction U(x, r) are highly anisotropic. Still at 
large length/time scales, U(x,t) becomes isotropic and 
logarithmic, which can also be confirmed by the asymp- 
totic expansion for large r: 



In' 



i 1 07 * 

lE + 2 In 2 - - In 5 - — \ + 0(l/r, 1/A 2 ). 



(15) 



Although the anisotropy in the short-range behavior of 
U may slightly alter the details of the RG flow, the qual- 
itative features at large length scales are expected unaf- 
fected. Thus it is concluded that for any value of finite A 
the system undergoes the superconductor-insulator tran- 
sition of the BKT universality class. 

Figure || displays the schematic phase diagram on the 
Kq-Ki plane. Although the precise phase boundary, i.e., 
the detailed behavior of Kq(X), in general depends on 
the microscopic length scales such as the (imaginary) 
time slice <5r, the universal (large-length scale) behav- 
ior should not be affected. Thus the phase boundary 
between the superconducting and the insulating phases 
is concluded to belong to the BKT universality class. As 
A (= Ki/Kq) is increased, the vortex creation energy also 
grows monotonically, resulting in the decrease of -Ko(A) 
toward the limiting value 2/ir. (Note that the value 2/ir 
corresponds to the zero vortex fugacity or the infinite 
vortex creation energy.) This lowering of the transition 
point and the resulting enhancement of superconductiv- 
ity reflects the role of the junction capacitance: It tends 
to suppress the quantum fluctuations induced by the self- 
capacitance, as manifested in the phase representation 
given by Eq. (|l4|). There has been few experimental 
studies of quantum phase transitions in ID Joscphson 
jucntion arrays. Recent interest in ID tunnel junction 
arrays |l4| ] and advances in submicron fabrication tech- 
niques strongly motivate experimenters in this field. 
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FIG. 2. Plot of the vortex interaction U(x,t) for (a) 
A 2 = 0.1 and (b) A 2 = 30. 




chain at zero temperature. As A gets increased, K§ decreases 
toward the limiting value 2/ir. 
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